Magnetic Excitations of Coupled Spin Ladders: a Quantum Monte Carlo Study 
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We calculate the magnetic excitation spectrum in the stripe phase of high-T c materials. The 
stripes are modeled as coupled spin- 1/2 ladders and the spin dynamics is extracted using Quantum 
Monte Carlo (QMC) simulations, which can capture the strong quantum fluctuations near quantum 
critical points of coupled spin ladders. We find a characteristic hourglass magnetic excitation spec- 
trum with high-energy peaks rotated by 45 degrees compared to the incommensurate low-energy 
peaks in good agreement with the experimental data. The excitations are investigated quantitatively 
as a function of interladder coupling, ladder width, and domain formation with stripe disorder. 
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Neutron scattering experiments are powerful probes 
of the bulk spin excitations of magnetic materials, and 
has been applied extensively to the cuprates since their 
discovery [lj. However, only recently have observations 
of universal fluctuations been reported0,0], providing a 
crucial new piece of information for the general under- 
standing of how these materials evolve from an antiferro- 
magnetic (AF) Mott insulator to a d-wave superconduc- 
tor as electrons are extracted from the CuC>2 planes. 

There is increasing evidence that competing interac- 
tions in the cuprate materials result in spatially inho- 
mogeneous electronic spin and charge ground states in 
a significant region of the phase diagram, and it is im- 
portant to determine whether experimental bulk probes 
are in agreement with such inhomogeneous solutions 
in general.and one-dimensional stripe configurations in 
particular While it might be tempting to model the 
magnetic response by checkerboard order 0, it is well- 
known that such textures cannot reproduce the dataQ. 
A recent neutron study of Lai.4sNdo.4Sro.i2Cu04 has 
confirmed that indeed the magnetic order is collinear and 
modulated one-dimensionally consistent with stripes 0. 

Recent extensive neutron scattering experiments on 
several cuprate materials have been performed to map 
out the details of the magnetic fluctuations over a 
large energy range 0, 0, 0, EJ- The compound 
Lai.875Bao.i25Cu04 (LBCO) is known to support static 
stripe spin and charge order below T = 50K as indicated 
by super lattice Bragg peaks in both sectors The low- 
energy magnetic excitation spectrum of this material 
(at T > T c ) consists of inwardly dispersing incommensu- 
rate spin branches. In L^-zSr^CuC^ (LSCO) the mo- 
mentum position of similar incommensurate modes are 
known to scale linearly with the doping up to ~ I a 
property that is naturally explained within the stripe 
scenario [l2l Hjjj . For LBCO the low-energy branches 
merge into the "resonance point" at (7r,7r) at an energy 
close to 50 meV. At even higher energies, constant en- 
ergy cuts reveal that the excitations disperse outward and 
form a square or ring feature with four prominent peaks 



rotated by 45 degrees compared to the low-energy incom- 
mensurate modes. The resulting spin response exhibits 
a characteristic hourglass shape, and is remarkably simi- 
lar to the data obtained in LSC OI14l and YBa2Cu306+ic 
(YBCO) H S ES E E El ET This supports the 
notion of a universal spin spectrum|l9l| modulo various 
quantitative differences, particularly at low energies re- 
sulting from e.g. material dependent exchange couplings, 
quasiparticle damping effects, and possible existence of a 
doping and temperature dependent spin gap. 

Theoretically, several groups have used various spin- 
only models to compute semiclassically the magnetic ex- 
citations in stripe states with magnetic long-range or- 
der (lro)E m nn. It has, however, been pointed 

out that strong quantum fluctuations are important for 
obtaining the correct high-energy response originating 
from the quasi-one-dimensional spin excitations of the 
ladder |2l 1221 l23l l24j. even though, in certain limits, the 
magnetic fluctuations can closely resemble the semiclas- 
sical results [i^. The importance of the charges in de- 
scribing the spin response remains controversial, but is 
expected to affect the details at low energy as mentioned 
above. The charge degrees of freedom have been in- 
cluded within phenomenological models [26j|. and within 
time-dependent Gutzwiller approximation of the Hub- 
bard model 24] . The coexistence phase of static stripes 
and cZ-wave superconductivity was studied in Ref. |53 us- 
ing a mean-field+RPA approximation. 

In this paper, motivated by the discovery of universal 
spin fluctuations, we calculate the dynamical structure 
factor 5(q, w) with a powerful QMC method able to treat 
both the ordered and quantum disordered regimes. We 
discuss effects of interladder coupling, ladder width, and 
disorder. As opposed to previous microscopic spin-only 
models, we do not rely on semiclassical linearized spin- 
wave ^lIlclEl Hll, or one-triplonp^ approximations. 

The spin dynamics is qualitatively different in isolated 
even and odd leg spin-1/2 ladders |28|: whereas even leg 
ladders display a spin gap and short-range spin-spin cor- 
relations, the odd leg ladders are qualitatively similar to 
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the spin-1/2 chain with gapless quasi-long-range correla- 
tions. For coupled ladders, the odd leg ladders support 
an ordered state for any interladder exchange coupling 
Jb, whereas coupled even leg ladders exhibit a quantum 
critical point (QCP) at some finite value of Jb 29] . For 
the cuprates, at a doping near x = g the charge is mod- 
ulated with a period close to four lattice constants, half 
the period of the associated spin modulations. Hence, 
the associated stripe phase may consist of coupled 2-, 3-, 
or 4- leg spin ladders containing 2, 1, and spin-empty 
sites, respectively. In the following, the spins on the Cu 
atoms in the CuC-2 planes are modeled by the Heisenberg 
Hamiltonian defined on a square lattice 

H = Y,JiiSi-Sh (1) 

(ij) 

where Si is the spin-1/2 operator at site i, and Jij de- 
notes the exchange coupling between sites i and j. Uni- 
directional stripes parallel to the y-axis are modeled by 
couplings Jf, across the hole-rich quasi-one-dimensional 
regions and J a elsewhere [3(| . The Hamiltonian and 
the assumed stripy configurations of exchange couplings 
is the resulting effective spin-only model after integrat- 
ing out the charge carriers. Whereas J a > 0, the sign of 
Jb depends on the period of the stripes: J& > for site- 
centered stripes and J& < for bond-centered stripes |2fl|. 
The choices accommodate the anti-phase property of the 
stripes believed to originate from a lowering of the ki- 
netic energy of the carriers. Below, we study L x L lat- 
tices with periodic boundary conditions, and use tetrag- 
onal units with wave vectors q x , q y parallel to the Cu-0 
bonds, and energies are given in units of J a . We typically 
set L — 64 allowing us to study momentum resolved cor- 
relations functions within QMC. 

The QMC simulations were performed using the 
Stochastic Series Expansion (SSE) method [3~J| with 
directed-loop updates |32|. In order to extract imaginary- 
time correlation functions the individual vertices in 
an SSE configuration were each assigned a random 
imaginary-time location in such a way that the as- 
signment did not alter the imaginary-time ordering se- 
quence of vertices[33j. The (Sf (t)S'J(O)) correlation 
function was read out using a mesh with 200 grid points 
along the imaginary-time axis and Fast Fourier trans- 
forms were used to obtain the momentum space points. 
The imaginary-time data was finall y c ontinued to real- 
frequency using the MaxEnt method [34j with a flat prior. 

Fig^shows superimposed results for S(c[, uj) as a func- 
tion of energy uj and momenta along cuts (q x ,i^) and 
(tt, q y ) for coupled 2- (a), 3- (b) and 4-leg ladders (c), all 
with | J(,| = 0.1 for easy comparison. For this value of Jf, 
the coupled 2-leg ladders are in the quantum disordered 
regime and displays a gap as seen in^a), whereas mag- 
netic LRO and Goldstone modes exist for the coupled 
3- and 4-leg ladders. Fig. [Ja-c) clearly show an hour- 
glass spin spectrum consisting of anisotropic low-energy 
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FIG. 1: (Color online) (a-c) S(q x , tt, uj) vs q x superimposed 
with S{ir,q y ,uj) vs q y showing the characteristic X-shaped 
hourglass spin excitation spectrum for coupled 2- (a), 3- (b), 
and 4-leg ladders (c) with \J b \ = 0.1 and J a f3 = 80. (d) 
S(ir,q y ,uj) vs q y from (c) showing the high-energy quasi ID 
dispersion along the 4-leg spin ladder (note different scale). 
In (b,c) the amplitude of peaks at (1 ± j)tt are reduced by an 
order of magnitude compared to the main peaks at (l ± j) tt. 

spin modes and a dispersive high-energy branch charac- 
teristic of quasi-one-dimensional quantum paramagnets, 
as shown separately in Fig. ^d). Important features in 
Fig. Hare the energy lj„ of the saddle point in the disper- 
sion at (tt,tt) (resonance point), and the minimum gap 
A at (l ± i) tt. Fig. 121a) shows the resonance energy 
for the 2-, 3-, and 4-leg ladders as a function of the inter- 
ladder coupling \Jb\- As seen, w T increases with \Jb\, and 
the even leg ladders display a gap at Jb = as expected. 
Note that for the 3-leg ladder oj^ does not vanish identi- 
cally at Jb = due to the finite temperature. This agrees 
with cr-model calculations which give ui^ = ttT+..., where 
ellipses denote higher order logarithmic corrections [35j. 

The intensity anisotropy ratio Rj between the inner 
(toward (tt.it)) and outer (away from (tt, tt)) branches 
of the low-energy modes is another important quantity 
since the outer branches are not yet observable in ex- 
periments. We define an anisotropy ratio by Ri = 
[S(q B + qi,7t) - S(q B - g/,7r)] / S(q B + qi , tt) where S(q) 
is the energy integrated structure factor, and qi = 
\qb ^ Qx\ with q B = ?>tt/A one of the Bragg points. 
Fig. EJb-d) show Rj vs qj for the coupled 2-, 3-, and 
4-leg ladders for various values of Jb- Cl early Rj de- 
creases with increasing interladder coupling |2ja], and for 
the 2-leg ladders Ri remains small (< 15%) in the vicin- 
ity of the Bragg point for all Jb- In fact, based only 
on Ri one would expect the 4-leg ladders to best model 
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FIG. 2: (Color online) (a) Resonance uj-n and mode gap A at 
37r/4 vs interladder coupling Jt- (b-d) The anisotropy ratio 
Ri vs qi — \qB—Qx \ for coupled 2- (b), 3- (c), and 4- leg ladders 
(d). All results are obtained with J a /3 = 40 and L — 64. 

the experiments!! 0, E3, EH H, EE EE E3, El For 

large enough Jf,, Rj eventually becomes negative and 
the outer branches dominate the intensity. We have also 
studied coupled 6- and 8-leg ladders, and found excita- 
tions qualitatively similar to Fig. ^b-c) but with even 
larger anisotropy ratios between the inner and outer low- 
energy spin branches. 

Fig. |21 shows constant energy cuts of S(q, w) for the 
coupled 4-leg ladders with the same parameters as in 
Fig. n^c-d). Similar to experiments, with increasing en- 
ergy the incommensurate low-energy peaks merge at the 
resonance point around ~ 0.3 and the high-energy re- 
sponse consists of a diamond structure dominated by four 
peaks rotated 45 degrees compared to the ones at low- 
energy. In Fig. 01 we show the momentum integrated 
structure factor S(ui) for the same parameters used in 
Fig. n Besides the ui — Bragg peak (for the ordered 
cases), S(ui) is dominated by two characteristic peaks 
corresponding to the saddle point in the dispersion at 
the resonance point, and a peak at the upper band edge 
of the dispersive branch along the ladders (uj ~ 2 J a ). 

Important open questions relates to the effects of disor- 
der and stripe domain formation on the magnetic excita- 
tions, and the origin of the spin gap. In LSCO magnetic 
LRO is observed at low doping x whereas 
at higher doping a spin gap phase, which may be related 
to superconductivity, opens up. YBCO and BSCCO ex- 
hibit a large spin gap, i.e. the stripes remain fluctuating 
at all doping levels possibly related to a smaller interlad- 
der coupling Jf, leaving these materials in the quantum 
disordered regime. Within a phenomcnological approach 
Vojta efc a].|3j| recently studied the spin susceptibility 
in the presence of slow dynamically fluctuating charge 
stripes. As a further complexity there is also evidence 
from muon spin relaxation experiments that LSCO and 
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FIG. 3: (Color online) Constant energy cuts of S(q,uj) for 
coupled 4-leg ladders; Jb = —0.1, L = 64, and J a /3 — 80. 
Similar results can be generated for coupled 2- and 3-leg lad- 
ders but with reduced anisotropy ratios Ri (Fig. |5Jl. 

YBCO exhibit a spin-glass phase with local 'static' short- 
range order. 

In the rest of this paper, we wish to study a simple 
disorder scenario which simulates both random stripe po- 
sitions and finite sized stripe domains. Domain forma- 
tion was an assumption when symmetrizing the plots in 
Figs. Hand ED and should be a natural consequence of 
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FIG. 4: (Color online) S(u) for coupled 2- (black) 3- (red) 
and 4-leg (blue) ladders with parameters similar to Fig. 
Dashed lines show the results for 400 x 12 sites with ordered 
(black) and disordered (red) 4-leg ladders (Jf, = —0.1). 
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random quenched disorder |38|. In the following we sim- 
ulate this kind of disorder by studying systems of size 
400 x L where L is an assumed average domain size, and 
the stripe positions are disordered with a flat distribu- 
tion in the interval [1,7] (mean 4). Clearly this approach 
is crude, but has been successfully applied to explain the 
quasiparticle spectral weight as detected by photoemis- 
sion spectroscopy^. In general, we find that this kind 
of disorder is very detrimental to the hour-glass excita- 
tions unless L is small enough to induce a spin gap. For 
example, for 400 x 12 systems with parameters similar 
to Fig. |3| we show in Fig. [5] the result for ordered (a) 
and disordered (b) stripes. The disordered case is aver- 
aged over 4 different configurations which is enough for 
these large systems. The corresponding S(u>) are shown 
in Fig. Fig. [3 demonstrates a case of domain-induced 
spin gap which is a property of the individual 4 x 12 clus- 
ters that (when coupled) constitute the system (a gap is 
negligible in the 64 x 64 system (Fig.EJ, and in 400 x 12 
systems with ~ J a (not shown)). Interestingly, as seen 
from Fig. OJb) , the stripe disordering not only smears the 
excitation spectrum resulting in diffuse 'legs of scatter- 
ing', but also strongly reduces the spin gap. In this case 
the resulting pile-up of weight right above the spin gap 
is clearly unrelated to (but may be enhanced by|2jj) su- 
perconductivity. We expect this scenario to be most rel- 
evant for LSCOQ, whereas in YBCO, where the inten- 
sity is peaked near (ir, j r), fluctuating charge stripes[3^ 
or superconductivity [27j may be important. For this ma- 
terial, however, it remains controversial whether homo- 
geneous models are better start ing poi nts, at least for the 
optimally-to overdoped regime |40j,|41|. 

In summary, we have calculated the spin excitations 
in the stripe phase within a Heisenberg model of coupled 
spin ladders using QMC, and investigated the effects of 
interladder coupling, ladder width and simple types of 
disorder. In agreement with experiments, we find low- 
energy spin-wave like excitations with intensity strongly 
dominated by the inner branches for low interladder cou- 
pling. These merge at the (jr, ir) point and rotate into 
the quantum excitations characteristic of the quasi-one- 
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FIG. 5: (Color online) "Boomerang" excitation spectrum for 
a 400 x 12 system of ordered (a) and disordered (b) stripes. 
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